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An effective action is obtained of a Bose gas in the bulk separated into two regions by a strong 
^ ' external potential depending on the single coordinate. The main attention is focused on the 

^ ' relaxation of the difference between phases of the weakly coupling condensates of the different 

^ bulk domains separated from each other by the external potential. 

^ . 

The experimental realization of the Bose-Einstein condensation at ultralow temperatures in atomic vapors [1] 
provides an example of the systems in which the approximation of the weakly non-ideal gas is well applicable. This 
• fact is connected with the small density of particles in the systems concerned. The realization of such systems gives the 
possibility of the experimental investigations of the macroscopic quantum phenomenons of different types. Recently, 
the manifestations of the macroscopic quantum phase and its behavior are studied actively [1-8], [11]. 

The investigation of the kinetic phenomena due to the relaxation of the order parameter is of a doubtless interest 
for the study of macroscopic quantum phenomena. In the present work we study the spatial relaxation of the phase 
of the order parameter in the inhomogeneous Bose gas with the weak coupling between different spatial regions due 
to a barrier tunneling. 

I The system concerned in the present work is a Bose gas of the small density in the bulk separated into two regions 
by a strong external potential U ('r') which depends on the x coordinate and does not depend on ±_ = (jj,z), 
I ' i.e., U (~f*) = U (x). For simplicity, we assume that the external potential has a rectangular shape U (x) — Uq for 
' O ■ —d < X < d and U (x) = beyond this region of the x coordinate. The height of the external potential Uq is supposed 
^ to be the largest energy parameter in the system, in particular, Uq » ^ where /i is the chemical potential = nA, 
A is the scattering amplitude of Bose particles, and n is the density of the Bose gas. Due to this assumption the 
interaction between particles of the Bose gas in the region of the influence of the external potential, i.e., in the region 
— d < X < d, can be neglected. 

The left and the right domains of the bulk are supposed to have the same temperature T smaller than the Bose- 
QQ . condensation temperature Tc- For simplicity, we suppose that in the initial state the density of the left and the right 
' domains of the bulk are the same, but the phases of the left and the right Bose condensates are different. This 
' assumption means that at the initial time moment there is nonzero current from one side of the bulk to another. This 
[ non-equilibrium initial state will relax to the equilibrium one which has the same densities and the same phases of 
I the condensates for both sides of the bulk. This relaxation process is studied in the present work. 
I Note that in the case of superconductors the similar initial non-equilibrium state results in the Josephson oscil- 
"Y"; ' lations with a small damping [9], [10]. The essential distinction of the Bose gas consisting of neutral atoms from 
' superconductors in which the Cooper pairs represent the charged objects is the presence of the gapless excitation 
^ . spectrum. The presence of low energy excitations results in the change of the character of the initial state relaxation 
I ' making the relaxation essentially faster. 

^ In the case of the homogeneous in y-z plane difference between the phases of the condensates the excitations 

Q ' radiated during the relaxation process have the one-dimensional character. The one-dimensional character of the 
Q radiated excitations results in a divergency of the small momentum correlators of these excitations. Hence, the 
' consideration of the relaxation process to second order in the tunneling amplitude is not correct in contrast to the 
case of superconductors [10] and requires more accurate analysis. 
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I. PERTURBATION THEORY IN THE TUNNELING AMPLITUDE 

As the introduction, we consider the perturbation theory in the amplitude of the tunneling under the external 
potential barrier. The Hamiltonian of the Bose gas with the tunneling term is 



H = J d^i 



(^) (-^^' - m) ^^(^) + ^ (V> ( r^) V^(^))'} + (1.1) 
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The form of the kernel of the integral operator Tyi-^yi-^ in the case of the external potential independent of can 
be taken as 

TV^T^^ = 72;i,2;2'5 (^±1 - ^-L2) where Xi > and X2 < (1.2) 

The kernel 7xi,x2 is localized near xi = X2 = within the coordinate x of about a where a < being the 
correlation length of the weakly interacting Bose gas (, = (nA) ' . In the case of the slow variation of the operator 
•tp (T*) compared with the scale ^ the kernel 7xi,x2 can be taken in the form 

7xi,x2 = 70^ (a;i) S {X2) (1.3) 

The constant 70 is the small quantity due to the large magnitude of the external potential and obeys the inequality 
7o << n^. The tunneling amplitude (1.3) will be obtained below and has the form 70 = ^n£, (>£'sinh {2xd)) ^ where 
K = \/2mU(j. Form (1.3) of the tunneling amplitude differs from that used in the work [11]. This difference is due 
to the behavior of the Fourier components of the tunneling amplitude 7fcij,,fe2x for the small momenta kix,k2x- In 
contrast to the present work, in the work [11] the tunneling amplitude 7fci^,/c2x vanishes "fki^M^ provided that the 
momenta kix,k2x tends to zero kix, k2x — > 0, this amplitude being nonzero for the condensate-condensate tunneling. 
This distinction has the far-reaching consequences, namely, nonzero value of Jki^^Mx for k\x-,k2x — * results in the 
divergency of lower orders of the perturbation theory in the tunneling amplitude. This problem is discussed in detail 
below. 

The characteristic frequency of the phase difference evolution between the right and left condensates is proportional 
to the small tunneling amplitude 70. The characteristic momenta of particles corresponding to this frequency are 
small compared with the scale of the order of 1/^. This means that the slow relaxation process of the phase difference 
is described by the slow variation in the spatial and temporal Fourier components of operators if) a {'r^,t). Index a 
is equal to a = 1, 2 where 1,2 denotes the right-hand (x>0) and the left-hand (x<0) sides of the bulk, respectively. 
Therefore, the model form of the kernel 7xi,x2 given by Eq. (1.3) corresponds to the purpose of the work, namely, to 
describe the relaxation of the phase difference. 

Due to the existence of the condensates in the left-hand and the right-hand sides of the bulk the operators V'a ("^, t) 
can be represented in the well-known form 

iia{^,t) = ^a{t)+Xa{^,t) (1.4) 



*o (t) = v^e^'^-^*) (1.5) 

The fields {t) are the condensate fields for the right-hand {a = 1) and the left-hand (a = 2) sides of the bulk. The 
value no is the Bose condensate density and (pa (t) are the phases of the right-hand (a = 1) and the left-hand {a = 2) 
Bose-condensates. The operators Xa {'r', t) are the annihilation operators of non-condensate particles. As has been 
written above, the main problem is the investigation of the time evolution of the magnitude A(fi (t) = ipi {t) — ip2 {t). 
For this purpose, we will obtain the effective action for the phases (pi {t) and (p2 (t). 

At first, we consider the perturbation theory in the tunneling amplitude for the small phase difference Aip (t). The 
diagrams for the renormalized tunnelling amplitude of the condensate particles have the form of the vertex with two 
external lines (as for the self-energy part) corresponding to the condensate field. The quantities corresponding to 
the incoming and outgoing dotted lines are ^'q, and '^a, respectively. To analyse the relaxation process, we use the 
Keldysh diagram technique with the Schwinger-Keldysh two-time contour [14], [15]. 

The reason of the impossibility to consider the lower orders of the perturbation theory in the tunneling amplitude 
is the divergency of the lower order diagrams when the external frequency uj tends to zero w — > 0. To demonstrate 
this, we consider the frequency Fourier transformation of one of second-order tunneling vertex 

rfi^=7o5(w)7o (1.6) 

where the function g (w) denotes the Green function of the non-condensate fields [13] with the momentum li ± = 
and the same x-coordinates 

g {to) = G\ iui, =0;xi = X2 = d) = G2 ( w, "^j. =0;xi = X2 = —dj (1.7) 
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The Green functions Ga {t — t',T^ — "f*') (a = 1, 2) are the Belyaev Green functions of the non-condensate particles 
of the homogeneous interacting Bose gas. The momentum k j_ corresponding to the Green function g (cj) vanishes 
k ± = due to the homogeneity of the condensate field in the y and z-directions. The Green functions Ga (^oj, k^ 

and, correspondingly, g (w) in the Keldysh diagrammatic technique has three components in triangle representation, 

namely, retarded G'^^\ advanced G*^'^^ and kinetic G^-^^ Green functions [14], [15]. In the case of small frequencies 
Ll) « n the retarded Green function g^^^ (w), for example, can be written in the form 

.^^^M-/%f . (1.8) 



{ui + i5) — c^p^ 

The small momentum asymptotic of G in Eq. (1.8) is used due to the smallness of the external frequency oj {oj « fx). 
The main contribution into integral (1.8) is provided by the momenta Px ~ w/c << /i/c. Note that integral (1.8) 
has a singularity and, therefore, the magnitude of this integral is determined by the bypass around the pole. The 
calculation of the integral (1.8) gives 

»™M = -.^^ (1-8) 

This behavior of the fimction g'-'^' (uj) means that for the sufficiently small magnitudes of lu, the consideration of some 
lower orders of the perturbation theory is not correct and the analysis of diagrams in all orders is necessary. 

The diagrams of odd order of this perturbation theory results in the renormalization of the bare tunneling amplitude 
7o. The sum of the odd order diagrams having one incoming and one outgoing line is denoted by the vertex r^p 
(a ^ (3) where, as above, the indices a, (3 = 1,2 denote the right-hand or the left-hand sides of the bulk. The sum 
of the odd order diagrams having two incoming external lines or two outgoing external lines is denoted by the vertex 
Ta/j [a 7^ /?) according to the notations of the self-energy part in the Belyaev diagram technique for the Bose systems 
[12], [13]. The term of the effective action corresponding to this renormalized tunneling amplitude has the form 



Stunn = j> j> dtdt'^^ (t) Tap {t, t') ^p it') + 

i dtdt' {t) Tap (t, t') ^p it') + (t) Tal3 (t, t') '^p it')} 



The time contour in this expression is the Schwinger-Keldysh two-time contour. Note, that in the even orders of the 
tunneling amplitude the renormalized vertex Tap has the coincident external indexes a = (3. These vertexes give the 
contribution to the self-energy part of the particles in one of the two subvolumes, in particular, the imaginary parts 
of these vertexes describe the relaxation process. 

The renormalized vertexes connected with the dressed tunneling amplitude can be expressed by means of the Green 
functions taking into account the tunneling 

Til = 70G2270 (1-10) 
Til = 70-F2270 

ri2 = 70 + 70G2170 

ri2 = 70-F2170 

The system of the Dyson equations for Gaji and Fap can be written as 

G22 = G2 + G270G12 + F270F12 (1.11) 

F22 = -F2 + -P270G12 + G2J0F12 

G12 = G170G22 + ^170-^22 
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Fi2 = G170-F22 + -F170G22 

Here Ga and Fa are the Bclyaev Green functions of non-condensate particles of the homogeneous interacting Bose 
gas. The Green function Ga has one incoming and one outgoing end and the Green function Fa has two outgoing ( 
or two incoming) ends. The Green function Ga denotes the Green function Ga with the opposite frequency oj and 
momentum i.e., Ga (w, = Ga (— w, — '^)- This system of equations is written in the assumption of the small 
magnitude of the difference between phases of the left-hand and right-hand condensates. In the case of the small 
momenta, which is the most interesting for us, the Green functions Ga and Fa obeys the equality Ga = —Fa- In the 
case of the small momenta the solutions of Eqs. (1.11) can be simplified and the tunneling vertexes take the form 



9 {^) 7o 
I-470VH 



i 11 - --i — ..^2„2 (. .^ (.1-12) 



Fii = 



ri2 = 70 + 



ri2 = 



9 (^) lo 

1 - 47252 {^) 



1-470V M 

VM7o' 



1 - 472^2 (^) 

In the limit w these vertexes can be written as 



" 45 H 2n ^ ' 

111 -. — = w 



1 . c 
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Ag (w) 2fi 



ri2 -2^0 
ri2^ - 2'^° 

Note that the bare tunneling amplitude 70 is eliminated from the renormalized vertexes Fn , Fn, and just these 
vertexes describe the relaxation of the phase (difference). This fact results in the large value of the relaxation rate. 

The more rigorous and detailed consideration of the problem in the case of an arbitrary phase difference for the 
infinitely large size of the system in the x-direction is represented below. The assumption about the size of the system 
is accepted to focus the attention on the consideration of the relaxation process. 

II. THE TUNNELING ACTION 

Below the approach to the problem of the phase difference relaxation without using the perturbation theory in 
the tunneling amplitude is developed. For this purpose the effective action for the boundary values of the Bose field 
phases is obtained by the integration over the volume components of these fields with the fixed values of them at 
the left-hand and the right-hand sides of the external potential boundary. Integrating over these fields we take into 
account the inhomogeneous behavior of the condensate fields in the left-hand and the right-hand sides of the volume. 
The necessity of the consideration of the inhomogeneous condensate behavior even for the small momentums of the 
excitations is conditioned by the fact that the self consistent inhomogeneous condensate additional field makes the 
situation in some sense similar to the case of the scattering on the external potential having the discrete energy level 
which is very close to zero. This can results, for example, in the transmission coefficient of the order of unit including 
the case of the very small compared with n energies [16]. The obtained effective action contains the relaxation part 
which is proportional to the first power of the frequency of the phase. For the small frequencies the relaxation part is 
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found to be much larger than the usual kinetic part of the phase dynamics which is proportional to the second power 
of the frequency. The consideration is developed in the framework of the Schwinger-Keldysh technique with the time 
reversing contour [14], [15]. 

The generating functional of the system can be written in the form 

Z = j DijDi)e''^^+^i'i (2.1) 

The action S of the inhomogeneous non-equilibrium Bose gas is given by 

S = jdtjd'r 1^ [idt +tJi-U {x)] V; - i (^V^) (^V) - \ (V^V')' I (2.2) 



and term Sj reads as 



Sj = j> dt J d\ + j^) (2.3) 

Here j and j arc the infinitely small sources. The Planck constant h and mass m of a Bose particle are set equal to 
unity h = m = 1. 

Due to the large magnitude of potential Uo in the region xe [— rf, rf] the interaction between Bose particles can be 

neglected and the integral for Z over the fields ip,tp takes the Gaussian form in this region of the x-coordinate. In 
this connection the integral for Z over the fields tp {x,l^ ±;t) ,tp {x,!^ ±;t) for xG [—d,d] with the fixed magnitudes 
of these fields at the boundary of the region of the nonzero external potential can be calculated. Moreover, the 
characteristic scale of the i/i-ficld variation in the time and in the ±_ space is supposed to be much larger than 
1/Uq and I/^/Uq, respectively. It is these fields that describe the slow relaxation process possessing the characteristic 
frequency proportional to the small tunneling coefficient. Due to this assumption the tunneling amplitude of the ip 
fields has the local character for the time variable t and space variables T^±. 

The magnitudes of the '(/'-fields at the boundary of the region of the nonzero external potential are denoted as 

V'l (x = d,T'±;t) = V.i i^±;t) (2.4) 
V'2 [x = -d, 't±;t) = 'ips2 (^±; 

V^i {x = d, T*_l; t) = ^p^^ {~tr, t) 
■02 {x = -d, T±; t) = ■0^2 (^-l; *) 

The calculation of the integral for Z over the ^/i-fields in the region xG [—d, d] gives 

Z = j £)V-DV'e'(^'""+^="'--f+'^-''' (2.5) 

The part of the action Syoi is 

S,oi = jdt j dx j d\^l^{idt + tJi)^-\ iy^) (^V) - ^ (V^V-) ' I (2.6) 

(7(x)=0 

where the region of the integration over the x-coordinate is the sum of the regions -L<x<-d and d<x<L. 
The part of the action Ssurf can be represented as a sum of two terms S^J^^^ and s'f^^p^ [17] 

Ssurf = 8^"^ + S^'^'T^ (2.7) 

^suTj surf ' surj V / 

sitf = -^coth{2^d)fdt J V.i(r^±;i) l' + l^.2(r^±;t) P} (2.8) 

^S"^ = 2sinh'(2xd) / / '^''■^ ■f^^^ *^ + *^ ^^-^^ 
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where the magnitude k is x = \/2Uo. Later on we assume that the width d of the potential barrier is sufficiency 
large so that the inequality >cd » 1 takes place and, thus, sinh {2>cd) ~ cosh (2>fd) w ^e^^'' >> 1. The term 5^^" ^"^ 
describes the tunneling between the right-hand and left-hand sides of the bulk. Prom the form of this tunneling part 
of the action it can easily be seen that the tunneling amplitude for the condensate particles and for the low energy 
non-condensate particles has the same nonzero magnitude in contrast with [11]. 

Our goal is the calculation of the functional integral for Z over the ■^-fields in the bulk, i.e., in the region xe 
(— L, — (i) U {d,L) with the fixed values at the boundary of the region of the nonzero external potential, and, thus, 
obtaining the effective action for the fields at the boundary tpsi, i>s2- 



III. THE SADDLE-POINT APPROXIMATION 

Later on, it is convenient to change the variables in the following way 

r ^ ^ r (3.1) 
1 

-t 
M 

where ^ = Thus, the system of units which we use later measures the quantities of the dimension of length in 

units of ^ and the quantities of the dimension of energy in units of /i. 

Moreover, it is convenient to represent the fields ^ in the modulus-phase form 

V' = pe**'; V' = pe-*'^ (3.2) 

The integral over the fields ip in the right-hand and the left-hand sides of the bulk can be calculated within the 
framework of the saddle-point approximation. The saddle-point magnitude tpci {x,l^±;t) of the field ip obeys the 
Gross-Pitaevskii equation with the boundary conditions 

tpci {x = d, T±;t) = tpsi (^_l; t) (3.3) 
V'cZ {x = -d, 'f'±;t) = ips2 {'r'±;t) 

In the present consideration we use the saddle-point approximation with respect to the modulus of the field tp, i.e., 
p. The fluctuations of the modulus field p in the vicinity of the saddle-point magnitude p^"^^ are taken into account 
within the framework of the perturbation theory. At the same time, the phase (p is considered exactly. The reason 
of such consideration is connected with the one-dimensional character of the excitations which are radiated during 
the relaxation of the difference between the phases of the left-hand and right-hand condensates. The saddle-point 
magnitudes p^'^'^ of the fields p have the form 

pf^ (x, r^i^ ; t) = V^Fi (x, r i ; t) {x > d) (3.4) 

P2^^ (x, = y/n^F2 (x, ~t±;t) (x < -d) 

where iiq is the condensate density in the bulk far from the region of the nonzero external potential and 

Fi(x,l^±;t) =ta,iih{x - d + Xi{T±;t)) (3.5) 
F2ix,T'r,t) =tanh{x + d + X2i'r'j_;t)) 

The fields Xi^2 {'r^±',t) are determined by the boundary magnitudes of the field p^'^') (x = ±d) 

psi {'r'±;t) = p[''^^ (x = d, ~t±;t) = y/n^ttanh (Xi {T'±;t)) (3.6) 

Ps2 {'r^±; t) = p2^^ {x = -d, r*_L; t) = v^tanh {X2 {'f'±;t)) 

Due to the large value of the external potential Uo (or >c >> 1) the essential contribution into the integral for Z 
is given by the fields with the small values of the modulus at the boundary | 'ipsi,2 i'r^±',t) |<< ^/no. This means 
that the fields X-1^2 can be considered as the fields with the small values | Xi^2 |<< 1 and, thus, these fields can be 
expressed in the simple form in terms of the boundary values of the fields Psi,2 ("^_l; t) 

^1,2 {'r'±;t) = ^Lp5i,2 {'r'±;t) 
\/no 
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IV. DENSITY FLUCTUATIONS 



For the consideration of the fluctuations of the modulus, the field p is represented as 

p = p(-')+^,5p (4.1) 

The boundary conditions are taken into account in the definition of the field p'^'^'^ and, therefore, the fiuctuations of 
the modulus y/n^Sp should obey the zero boundary conditions 

5pi [x = d,T'±;t) = Sp2 {x = -d, r"^; = (4.2) 

Within the framework of the saddle-point approximation the inequality | 5p \« p'^^'V%A^ = ^ should hold for. 
Substituting (4.1) into the action (2.6) and taking into account the inequality | 5p \« /o'-'^'Vv^ ~ obtain 

Svoi = + SSyoi (4.3) 

where 

Sil\=noefdt J dxJd'r^S^-FHd,^)+F'-lF'-l(VFy^ (4.4) 



where the function F is obtained by Eq. (3.5) and the operator Hi reads 

1 



Hi = --Vl-l + ZF^ (4.6) 
The spectrum of the operator Hi can be found [19] and has the form 

Ao = 0, Ai = ^ , Afe = 2 + (4.7) 

Note that the zero mode sclf-cncrgy Aq of operator Hi should be eliminated due to zero boimdary conditions (4.2). At 
the same time the second discrete self-energy Ai and the self-energies of the continuous spectrum Afe are not changed 
essentially by boundary conditions (4.2). In Eqs. (4.4), (4.5) the fluctuations of modulus 5p are taken into account 
within an accuracy to second order due to inequality Sp << F and the phase is taken into account exactly. The; 
possibility of the consideration of the fiuctuations 6p within an accuracy to second order results from the gap character 
of the spectrum of these fiuctuations, in contrast to the spectrum of the fluctuations of the phase which is gapless. 

The contribution of the dcnsify fluctuations into the generation functional Z is determined by action SS^ol and has 
the form of the Gaussian integral over 6 p. The integration over Sp in Z can be performed and we obtain 

Z = l DXD^ei'-'+''''+'-^+'^) (4.8) 



dt 



d^l^d^ 



r" < 



{F (dt^))^ [Hi 0^) (F (dt^))^^ 



(4.9) 



The evolution of the difference between the phases of the condensates in the right-hand and the left-hand sides of 
the bulk is a slow process. The inverse characteristic time of this process is proportional to the tunneling amplitude 
7 = ^no (xsinh (2>fd))~^ and has the small value. Therefore, the momenta A; _l of the fluctuations corresponding to 
these fluctuation frequencies should have small value fcj,^ << 1 (or in our system of units where ^ = 1, fcj, << 1). As 
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a result, due to the gap character of the self-energies of operator Hi in the case of zero boundary conditions for 5p in 
the direction 'r^± we can neglect the term compared with Hi in the first summand in the action (4.9). 



The operator Hi can be represented in the factorized form 



Hi = -L+Li 



(4.10) 



where the operators Li, have the form 



Li = V, + 2F 



(4.11) 



Due to the possibility to represent Hi in the factorized form and to neglect in the first summand of Eq. (4.9), 
the action S'^'^'> can be written as 



S^'P^ =noe jdt J dx J (fr^h (1+) '{Fdtip) -ii^2(^<^) 



(4.12) 



!7(x)=0 



After the integration by parts over x in the second term of the action S^'^^ (4-12) this action can be transformed to 
the sum of two terms 



where 



Siol = "0<e^ fdt J dxj d^T"^ 
U{x)=0 



[Fif) H2 (Fif) 



(4.13) 



(4.14) 



Si:lf = -lnoefdt J rf^r^{[FVi(V.^i)] |^=^ + [FV2(V.<P2)] llZZi} 



(4.15) 



Ho 



-^Vl-l + F' 



(4.16) 



The term s'^'^lj of the action S^'f^ is a surface one. The region of the integration in this surface term is the surface of 
the right-hand and the left-hand sides of the bulk. 

Note that the operator H2 can be represented in the factorized form as well as operator Hi 



H2 = 2-^2-^2 



where 



L2 = -v^ + 
L+ = + 



F 

(VxF) 



(4.17) 



(4.18) 



It can easily be obtained that the operators Li and L2 obey the following equations 

LiH2Lt = ^LiL+L2L+ = ^V, (V^ - 4) V, 



(4.19) 



LiL+ = 2H2+4 



(4.20) 
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V. THE SURFACE TERMS OF THE ACTION 

The part of the action S^Ji which consists of the modulus terms alone and does not have the phase fluctuations 

can be transformed to the surface form. This can easily be seen by the substitution of equation (3.5) for F into the 
action S'^^j (4.4). Taking into account the smallness of X (X<<1), we obtain 

noefdt J da:ydV|F2-iF4-l(V.F)'-i(^^F)'| (5.1) 

t/(x)=0 

Note that taking into account the smallness of X and assuming the slowness of X compared with the scale ^ in the 
1^ ±_ space (the reason of this assumption has been discussed above), wc can neglect the term . 

The term of the action 5^^] depending on the fluctuation of the phase can be transformed to the surface form in 
the assumption of the slowness of the field compared with the scale ^. 

no^^jdt j dx j d^r:^{F^{dtip)) = (5.2) 



= nof jdt j dxj dW { {F^ - 1) (dtif)} = 

U{x)=0 

= no^^ jdt j d^ri_ {X {dt^,)} 



As a result, the sum of 5^°] and Ssurf is given by 

S^2+Ssurf=noe fdt I d\^ |x - - X [d^^,) + ^.^^^^^^^ cos (Ay,) | (5.3) 

where At^, = ^si ~ H^si is the difference between the magnitudes of the phase fields at the boundary of the region of 
the nonzero external potential. 

Due to the large value of n^^^ >> 1 the saddle-point approximation for the integral over X in the generation 
functional can be used. As can readily be seen, the saddle-point value of X is equal to Xq 

= - (5.4) 

H 

After integrating Z over the fluctuations of X near the saddle-point value Xq and taking into account (4), we obtain 
the surface part of the action in the form 

C = fdtj d^r. [-1 iO,^.f + ^3inh(2xd) ^^^^^l ^'''^ 

Combining two surface terms iS^^^j (5.15) and S^^^j: (5), we obtain the surface term in the form 



VI. THE EFFECTIVE ACTION FOR THE PHASE 

Due to the smallness of the fluctuations of X we replace X with the saddle point value Xq in the fleld F and denote 
this field as F(°) where 
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F(o)(.)=tanh(^i|^+Xo) 



(6.1) 



The function F*^°^ is substituted in the action S^^J defined by Eq. (4.13) instead of F. The same procedure refers also 
to the operators L^, and Li. After this replacement the bulk part of the action S^^) can be written as 



j> dt j dx j d' 



C/(x)=0 



"(£(o)+)-\,.(o)a^^ 
(i^(o)^)ijf (FW^)-i(fW)'(^^^)\ 



(6.2) 



where the operators hf'"'^ and i?2°^ have the form of and i?2 with the substitution F F^^\ Since the action 

introduce the field 



'^voi ^^^^ form, that it is convenient to introduce the field r] instead of the field ip 



The transformation yields 



si:} = i> dt 



dx I d^r± 



Using (4.19) and (4.20), the action 5^^/ can be written as 



U(x)=0 
(V) 



-1 (V.r,) (mr^-) 



77- 



(6.3) 



(6.4) 



Si:j=2noei^dt 



dx / d'^r± 



{dtvf - h [V. (Vl - 4) Vj 7?- 



u{x)=a 



) (^-) 



(6.5) 



In the case of the slow spatial variation of the field i] along the x axis the action (6.5) can be transformed to 
S^:^l = ^noe / di J dx j Sr^ I {d,rif + (V^) n - \ (^^r?) (l + (^^r;) 

!7(x)=0 

The transition to the Fourier transformation in the plane gives 



(6.6) 



(6.7) 



The surface terms in the action S'fj^^^ (4-15) can easily be expressed in terms of the field rj. Using the connection 
between the fields and r] (6.3) having the form 



we express in terms of rj 



1 + (F(0))'\ 



(6.8) 



(6.9) 



The characteristic frequency t^o of the phase time evolutions is proportional to the tunneling amplitude loo ~ 7, 
where the tunneling amplitude 7 = (xsinh (2>cd))~^ . The tunneling amplitude 7 is the small value being much smaller 
than = F^^^ (±d) << 1, i.e., ji in the usual system of units. Thus, the characteristic frequency of the field rj is 
small and, therefore, Vj;?] is small compared with r] due to the connection w ~ fcx on the mass surface for the small 
uj (Eq. 6.7). In this case Eq. (6.9) goes over into 
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l+(F(o))2 



(6.10) 



Prom Eqs. (6.8) and (6.10) due to the inequality wq ~ ~ 7 << F^^^ (±c^) (^x^^ is the characteristic momentum 
corresponding to uq) we obtain 



<fi{d) = 2rj{d) 



(6.11) 



(p(L) =2r?(L) 



Vc,<^ (i) = 2Vxr? (i) 

Prom these equahties the surface terms in the Eq. (4.15) for 5^^^^ can be written in terms of the 7] field in the form 



(f(0))%V.<p 



{d) = [vV.v] {d) 



(6.12) 



VII. THE INTEGRAL OVER THE PHASE 

Using expressions (6.5-6.7) for action S^^^j, the integral over the field ip {x,l^ ±;t) with the fixed values (ps (l^^;<:) 
on the surface of the region of the nonzero external potential can be calculated. Note that the phases of the left- 
hand and the right-hand condensates are homogeneous in the 'r'± plane and do not depend on "f*_L due to the 
translational invariance of the system in these directions. Therefore, the surface phase fields (ps (l^x; t) depending on 
"f*_L correspond to non-condensate particles and can be considered as perturbations. We analyse the relaxation of a 
difference in phases in zero order of these perturbations and do not consider them in this work, i.e., we put (psi {t) 
and (ps2 {t) depending on the time variable alone. 

The action S^^^l (6-5) has the Gaussian form. In this connection, in order to calculate the integral over the bulk 
components of the field r] for the small oj and k , we should find the solution of the equation 

V = (7.1) 

as a function of the boundary values of the field rj and put this solution into the action S^j'^i . We supposed that the 
boundary values of the phase field (p, i.e., ipg and ip^, are known. In the case of the boundary magnitudes of the field 
(p and, hence, of the field rj independent on we look for the solution of this equation for k ± = 0, i.e., 

[u;'+Vl]v = (7.2) 

As was mentioned above, to describe the relaxation, we use the Keldysh-Schwinger technique with the two-time 
contour [14], [15]. In this technique the field r] has two components corresponding to two branches of the time 
contour (quadrangular representation). It is convenient to introduce the components of the field 77 as a sum and a 
difference of the two components corresponding to the upper and to the lower branches of the time contour (triangular 
representation) [14], [15]. In Eqs. (7.1), (7.2) the frequency w is a complex quantity u ^ iv + iT. The imaginary part 
r is the relaxation frequency of the excitations in the homogeneous Bosc gas. The value F depends on uj and for the 
small Lo this value is much smaller than lu, i.e., lo >> T [20]. Thus, the equation for the first component of the field 
T] has the frequency w in the form of w — > w* = w — and the equation for the second one has the frequency in the 
form u! ^ LO = UJ + ir. ^ 

The Keldysh Green function of this equation is denoted as D {x,0). We use the "triangular" representation. The 
solution of Eq. (7.2) with the fixed boundary conditions is 



^2+v|-tl + itl(v.F(°)) 
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V 



(ar)=5(a;,0)(5(0,0)) \(0) (7.3) 



the coordinates x=±d of the boundary we denote as 0. The Green function D {x, 0) in the "triangular" representation 
reads 

Da 



where Dn and Da are the retarded and advanced Green func;tions respectively, and Dk is the kinetic Green function. 
For the equilibrium state, the Green function Dk is connected with the functions Dji and Da by the equation 

Dk = coth (^) {Dr - Da) (7.5) 

The Green functions Dji and Da have the form 

^k(x,0) = /|^^^— — ^= (7.6) 
J 2n pl-(aj + iT) 

^ ^ 2 (w + iP) ^ ' 2{u; + ir) 



Da {x,0) 



dpx 



2n pI-{w- iV) 



2 



^ ' 2{uj + iT) ^ ' 2{io + iT) 



The inverse Green function 1? ^ in the "triangular" form is 

5-1 = -^^^3'^^' ) = ^^^^^_t - ) (7.7) 

The space derivatives of these Green functions give (x>0) 



V,Cfl(a;,0) =-ie*^('^+^r) (7.E 
y,DA{x,Q) =-ie--(-'r) 



Using Eq. (7.3), we can represent the surface terms (5.12) and (4.15) of the action S^^^^ as 

[v^xV] (0) = (7.9) 



VxD A (0) 
V,£>fl(0) V,£>k(0) 





^0 












2 



-2iwcoth(^) —iijo 
iw 



Thus, we have 



(7.10) 



(5(0,0))-^ = (^coth(-) (z),^(o,^o)^-z?-(o,o)) D-,\,,,)y (^^^^ 



[.V..] (0) = ( (0) (0) ) ( ^f^ (#) ) ( I ) (7.12) 
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Substituting (7.12) into (4.15), we obtain the following expression for the action S^'^lj^ in the "triangle" representation 

Si:lf = -InoeS^ Jdti^i 4>s, ) (5(0,0))"' ( ^^Jl (7.13) 

In the time representation the frequency oj is replaced by the time derivative idt- The fields $ and (j) are defined by 
the equalities 

^ = \{Vi+V2) (7.14) 

= ^ ((Pl - <P2) 

where (pi and </?2 are the phase fields determined on the upper and the lower brunches of the time contour, respectively. 

As a result, the action -S^^^j as an effective action Sgff which describes the relaxation of the difference between 
the phases at the boundary of the region of the nonzero external potential can be written in the "triangular" form as 



where 




(7.15) 



(7.16) 



Vi^s) = . .^^ cos(y,) (7.17) 

xsmh [2Kd) 

The index a = 1, 2 in Eq. (7.15) corresponds to the right-hand or the left-hand side of the bulk and A$s = $si — $s2, 
= (f)s\ — (f)s2- The value S± is the area of the section of the bulk along the y-z plane. 



VIII. THE DESCRIPTION OF THE RELAXATION 

The effective action (7.15) describes the relaxation of the difference between the phases at the right and the left 
boundaries of the region of the nonzero external potential. The generation functional Z can be written in the form 



(8.1) 



To simplify the consideration, we suppose that the temperature T is much larger than the characteristic frequency 

~ 7 = (>fsinh (2>i:d))~"'^. In this assumption, firstly, the term coth(^) in (_D(0,0)) (7.11) can be written 

as Secondly, the component of the phase field <j) describing the quantum fiuctuations is much smaller than the 
component $ {(f) » $). This approximation is conditioned by the large value of the parameter no^^S± » 1 (the 
system of units is such that the scale of the length is defined by the value ^). Due to the smallness of the field (j) 
compared with the field $ it is possible to expand the nonlinear terms V {A^g ± A^g) in a series in A(j)s to first order. 
Thus, we obtain 

r { (t>ea[^Oj'^ + {'^]^sa + '^sa[^^'^ - li^](t>sc.+i(l)soT(t>sa-\ 

Seff =no^-^S± / dt{ «=i,2 } (8.2) 

The term V (A$,,) is the functional derivative of the functional V [A$,,], i.e., V [A^g] = 6V [A$s] /S (A^^). 
We introduce the half-sum and the half-difference of the fields (j)sa and 
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* = ^ + $«2) ; <t>=^{ct>si+ct>s2) (8.3) 
(A$) = i ($,1 - $,2) ; (A0) = i (0,1 - </.,2) 

After the integration of the generation potential Z over the fields $ and ^ we obtain 

Z = j D {A(j3) D (A$) e^(ASe//+Si) (8 4) 

where 

z^^e// not.5iyo(t| +i(A0)T(A0)-2(A0)r [A$] J ^^'^^ 

The fields A$ and A0 are represented as a sum of the parts changing slowly and fast in time. We denote the fast 
parts of A<1> and A</) as (5$ and 50 and the slow parts as $ and 0. Thus 

A$ $ + (5* (8.6) 

A0 + (50 

The criterion of separating the fields into the fast and slow ones is the magnitude of the frequency. If the frequency 
of the field w > Wc where ujc = the field is fast, correspondingly, if the frequency w << Wc, the field is slow. This 
gives 

_ nislow) aifast) /n 7^ 

^eff-^eff +^eff 1^-') 

5(/-t)_„.35, /rfJ (<50)[2b^2 + iu;-r'[<I>]](<5a>)+ 1 



The integration over 50 yields 



where the operator L is 



4-^"'*^ = y dt [l (5$)] 'I (8.10) 



L=^df + \dt + V"m (8.11) 

Note that for the slow frequencies w << 1, the term can be neglected compared with the term jdt- 

For the generation potential, we obtain 

Z = ZslovjZfast (8.12) 

where 

Zsiou, = J D0£)$e'(^="""+^^-''°"") (8.13) 

Zfast = J D(5cl>)e^''* +^^- ) (8.14) 
Due to the form of the action S^"'^*^ (8.10) the generation potential Zfast has the Gaussian form and can be calculated 
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Zfast = Cideti 



(8.15) 



where the constant C has the form 



and can be included into a measure of the functional integral Zsiow The index 1 in deti denotes the product of the 



Fourier components corresponding to the fast frequencies. The value deti 
way 



deti 



= C2 exp {-Sp'Ln [(of) V" ($)] } 



can be calculated by the following 

(8.16) 



The symbol Sp" indicates the sum over the fast frequencies and the symbol Sp' indicates the sum over the slow 
frequencies. The inessential c-number constant C2 is equal to 



C2 = exp { -Sp"Ln 



and can be included into the measure of the Zsiow The value Sp'Ln ^G^'^ V" ($) 



is calculated as 



5p'in [(eg)) y" ($)]=!/"($ (t)) j ^_-L_^ = -4^(0)y"($) = -2F"($) 
Hence, the expression for the generation potential has the form 



-'fast 



Cideti 



exp <^ 2 / dtV" ($) 



This equation gives 



Z = Z^iowZfaet = / D(j)D^e \ -St ^ ^ J ^ >) 

The action 5^^^**^ has form (8.9). Integrating Z (8.18) over the field (j), we obtain 



(8.17) 



5.18) 



(8.19) 



1 



1 



(9*$) - V ($) 



- 2V" ($) 



Omitting the term (9t$) V ($) which can be represented as a total derivative, we obtain 

sifL = dt\^^^{dt<^f + {V - 2TV" ($) 



(8.20) 



(8.21) 



Note that the generation potential Z with action S^fJ^ (8.21) has the form corresponding to the quantum mechanics 
with the " imaginary time" in which the role of the " imaginary time;" plays the; usual real time. This action corresponds 
to the action which is obtained for the classical random process describing by the Langeven equation with the white 
noise [21] and having the form of the supersymmentic quantum mechanics. 

The Hamiltonian corresponding to the action iS'^*^^^ (8-21) is 
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t2 

H = -4T^^ + {V i^)f - 2TV" ($) (8.22) 
and the Schredinger equation corresponding to this Hamiltonian can be written as 

-5t* = (8.23) 
The Hamiltonian H can be represented in the factorized form 



^ = -(2T)^(^-^)(^+^) (8-24) 



d .\ f d 



where 

2T 



1.25) 



The Hamiltonian H has the minimum eigenvalue equal to zero. The eigenwave function corresponding to the zero 
eigenvalue is 

*o = exp ($)| =cxp|-^5^| (8.26) 

This can easily be seen from the form of the Hamiltonian (8.24). Due to the form of the Hamiltonian (8.24) the 
Schredinger Eq. (8.23) can be modified into the form of the Fokker-Planck equation [18] by the transformation 



^ ^expl J d^'A ($') I * 



(8.27) 



In this case the eigenwave function (8.26) corresponds to the equilibrium distribution function for the field $ which 
is formed as a result of the relaxation process. 

This research was supported by the Russian Foundation for Basic Research and by the Netherlands Organization 
for Scientific Research (NWO). 
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